Search for Z^ x and Zf 2 strangeonium-like structures 
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Theoretically, it has been presumed from an effective Lagrangian calculation that there could exist 
two charged strangeonium-like molecular states Zf x and Z* 2 , with KK* and K*K* configurations 
respectively. In the framework of QCD sum rules, we predict that masses of Zf x (KK*) and Zf 2 
(K* K*) are 1.85±0.14 GeV and 2.02±0.15 GeV respectively, which are both above their respective 
two meson thresholds. We suggest to put in practice the search for these two charged strangeonium- 

£N) , like structures in future experiments. 

t-H ■ 

O PACS numbers: 11.55.Hx, 12.38.Lg, 12.39.Mk 

<n: 

C3 . I. INTRODUCTION 

S: 

1 Recently, Liu et al. study the (/>(1020)7r + invariant mass spectrum distribution of Y(2175) — > 
</>(1020)7r + 7r~ and indicate that there could exist two charged molecular states and Z^ 2 , whose con- 
_ _ , figurations are KK* and K*K* respectively [lj. The molecular state is well and truly not a new concept 
but with a history, which was put forward long ago in Ref. Q and has also been predicted that molecular 
states have a rich spectroscopy in Ref. j^j. Since there is not any restriction for the number of quarks 
^4 inside a hadron, QCD does not exclude the existence of multi-quark states such as molecular states. In 
fact, some of the so-called X, Y, and Z resonances have already been ranked aspossible charmonium-likc 
molecular candidates. For example, X(4350) is interpreted as a D*D®* state [J, |5j; F(3930) is proposed 
(Nl ■ to be a D*D* state 0-0); Y(4140) is deciphered as a D*D* stateE^; F(4260) could be a XcP° or an 



X 



ujXd state [llj; Y74274) is investigated as a A,D s0 (2317) state [12j; Z+(4430) is suggested to be a D*D 1 
molecular state [13J. For more molecular candidates, one can also see some other reviews, e.g. Ref. 14 1. 
The two Z+ and Z+ resonances may shed light on studying strangeonium-like molecular states. Their 
■ properties like masses are important and helpful for searching them in future experiments. Unfortunately, 
quarks are confined inside hadrons in the real world, and the strong interaction dynamics of KK* and 
K*K* systems are governed by nonperturbative QCD effect completely. The quantitative calculations of 
hadronic properties run into arduous difficulties. However, one can apply the QCD sum rule method fl5| 
. (for reviews see and references therein), which is a nonperturbative formulation firmly based on 
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QCD basic theory and has been successfully employed to research some light four-quark states 
this work, we are devoted to predicting masses of Z+ and Z+ from QCD sum rules. 

The rest of the paper is organized as three parts. We discuss QCD sum rules for molecular states in Sec. 
|TT]with the similar procedure as our previous works 3, 27 1, where the phenomenological representation 



and the operator product expansion (OPE) contribution up to dimension ten operators for the two-point 
correlator are derived. The numerical analysis is made in Sec. Mil and masses of Z^ (KK*) and Z~^ 2 
(K*K*) are extracted out. The Sec. IIVI includes a brief summary and outlook. 



II. QCD SUM RULES FOR Z+ AND Z+ 2 MOLECULAR STATES 

An elementary step of the QCD sum rule method is the choice of interpolating current. Following the 
standard scheme [28(, strange mesons with J p = 0~ and 1~ are named as K and K* . In full QCD, 
interpolating currents for these mesons can be found in Ref. . One could construct the molecular state 
current from meson-meson type of fields. Thus, the following form of current could be constructed for 
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KK*, 

3kk* = (scn59c)(g' c '7' t Sc') I (!) 

where q and q' denote light quarks u and d, c and c' are color indices, and the quantum num- 
ber for the current is 1+. Lorentz covariance implies that the two-point correlator U^(q 2 ) = 
i J d 4 xe tqx (0\T[j IJ -(x)j I/+ (0)]|0) can be generally parameterized as 

IP"V) = {^f -^n«( 9 2 ) + t*L U W(q% (2) 

where H^\q 2 ) is pure vector and ±l(°)(g 2 ) is related to the scalar current correlation function. In phe- 
nomenology, the calculation proceeds by inserting intermediate states for KK*. Parameterizing the cou- 
pling of the state KK* to the current j^^* m terms °f tnc coupling constant A^ as (0\j^^,\KK*) = 
A^e^, the phenomenological side of IP"^ 2 ) can be expressed as 

where M K ^» denotes the mass of the KK* resonance, and so is the threshold parameter. The Lorentz 
structure g^ v gets contributions only from the spin 1 state, which is chosen to extract the mass sum rule. 
In the OPE side, n (1) (g 2 ) can be written as 

U^(q 2 ) = ds? [ f. + U^ cond (q 2 ), (4) 

J 4m 2 s — q 

where the spectral density is (s). After equating the two sides, assuming quark- 

hadron duality, and making a Borcl transform, the sum rule can be written as 

[\W] 2 e- M iK*/ M2 = f S ° d.sp 0PE e-^ M2 + £n« cond . (5) 

J Ami 

where M 2 indicates Borel parameter. To eliminate the hadronic coupling constant A' 1 ', one reckons the 
ratio of derivative of the sum rule to itself, and then yields 



Jim* tthjp) J / [ J 



The current for K*K* could be constructed as 

3k*K' = {sc^<lc)(q'c^ s c'), (7) 

with the quantum number + . Phenomenologically, the correlator II(<7 2 ) = if d i xe lqx (0\T[j(x)j + (0)] |0) 
can be expressed as 



m k*r* -y 2 71 J so s - 1 2 



n(g 2 ) = nir2 K * K " + - / ds 5--^ + subtractions, (8) 



where M K ,^» denotes the mass of the K*K* resonance, and \k*k* gi vcs the coupling of the current to 
the hadron (0\j\K*K*) = \k*k*- I n the OPE side, the correlator can be written as 

U(q 2 ) = r ds P PE( 2 S) + n cond (g 2 ), (9) 

J Ami s - q 2 
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where the spectral density is p^PE^ _ li m lfiPE^ Then, the sum rule can be written as 

p8Q 

\ 2 K , R ,e- M K***' M2 = / rf S/) ° PE e - s / M2 +M cond (10) 
Eliminating the hadronic coupling constant A^-*^*, one yields 

For the OPE calculations, we work at the leading order in a s and consider condensates up to dimension 
ten, utilizing the light-quark propagator in the coordinate-space 

- J^siSQo- ■ Gq) + y—^m q (gq(T ■ Gq)^ - ^ w ^(qq)(g 2 G 2 ). 

The s quark is dealt as a light one and the diagrams are considered up to the order m s . The spectral density 
can be written as p 0PE (.s) = pP crt (s) + p<™>(«) + p< s ~ 5 >(s) + p<«> <">(«) + pbfr-Gq) ( s ) + p {gi<r-Ga)^ + 

P<3 2g2 )( S ) +p(99><9 2 G 2 )( s ) +p <gg><flS<7-G S >( s ) + p {S S ){gqvGq) where ^pert^ p (gg> j p (qq)(ss) ^ p (gqo-Gq) ^ and 

p(s G ) are perturbative, quark condensate, four-quark condensate, mixed condensate, and two-gluon 
condensate spectral densities, respectively. They are 



3-2 15 tt 6 ' r w 2 10 tt 4 s ' ' w 2 10 tt 4 6 ' r w 3-2 5 ^ 2 



<5 9 ><9 2 G 2 ) f ^ _ (gg)(ff 2 G 2 ) ( ft ) (sfo .G,) f , _ 3(w)^'Ga) ( Sa >< OT - < r.G,> fs1 _ 3(ss)(gqa- Gq) 



6 n (i)cond = rn s (ss) 2 (qq) _ m s (qq} 2 (ss} (gqa ■ Gq)(gsa ■ Gs) (qq)(ss)(g 2 G 2 ) 

3-2 3 3-2 2 8 tt 2 3 2 -2 7 tt 2 ' 

for KK*, and 

^«=5^M- 4 > P™(s) = -H$m s s\ P ^(s) = ^m s s 2 , ,<*><*>(,)= <^ fl) 

n {qq){gia.G,) ( .) _ (qq)(9s<T-Gs) { ss)( g qa-Gq) ( _x _ (ss)(gqa-Gq) 

n n cond = m s (sg) 2 (gg) _ 2m s (qq) 2 (ss) (gqa ■ Gq}(gscr ■ Gs) (qq)(ss)(g 2 G 2 ) 

3-2 3 2 6 tt 2 3 2 .2V ' 

for 
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III. NUMERICAL ANALYSIS AND DISCUSSIONS 

Numerically, sum rules (j6|) and are analyzed in this section. The input values are taken as m s = 
0.10±g;g| GeV Q, (qq) = -(0.23 ± 0.03) 3 GeV 3 , (gqa ■ Gq) = m% (qq), (ss) = -(0.8 ± 0.1) x (0.23 ± 
0.03) 3 GeV 3 , (gsa ■ Gs) = m 2 , (ss), m\ = 0.8 ± 0.1 GeV 2 , and (g 2 G 2 ) = 0.88 GeV 4 0. Complying with 
the standard criterion of sum rule analysis, the threshold sq and Borel parameter M 2 are varied to hnd the 
optimal stability window. In the QCD sum rule approach, one can analyse the convergence in the OPE side 
and the pole contribution dominance in the phenomenological side to determine the allowed Borel window. 
Meanwhile, the threshold parameter ^/sq is not completely arbitrary but characterizes the beginning of the 
continuum state, and the energy gap between the groundstate and the first excitation is around 0.5 GeV 
in many cases of light mesons and nucleons. In a word, it is expected that QCD sum rule's two sides 
have a good overlap in the work window and information on the resonance can be reliably obtained. For 
instance, the comparison between pole and continuum contributions from sum rule (1101) for K*K* for 
y/so = 2.4 GeV is shown in the left panel of FIG. 1, and its OPE convergence is shown in the right panel 
by comparing the perturbative, two-quark condensate, four-quark condensate, mixed condensate, two- 
quark multiply two-gluon condensate, two-quark multiply mixed condensate, six-quark condensate, mixed 
multiply mixed condensate, and four-quark multiply two-gluon condensate contributions. Note that the 
perturbative contribution is almost as large as the (qq)(ss) contribution at M 2 = 1.5 GeV 2 (the ratio of 
(qq)(ss) to perturbative is approximate to 96%). Even if we choose some weak convergence criteria, e.g. 
the perturbative contribution should be 20% bigger than the second most important condensate, there 
is no standard OPE convergence at least up to M 2 > 1.8 GeV 2 (the ratio of (qq)(ss) to perturbative is 
approximate to 79% at M 2 = 1.8 GeV 2 ). On the other hand, the relative pole contribution is approximate 
to 53% at M 2 = 1.3 GeV 2 and descends along with the M 2 . The consequence is that it is not possible to 
find a region where both the OPE normally converges and the pole dominates over the continuum. The 
problem with the sum rule is that the perturbative contribution is smaller than the four-quark condensate 
contribution while the pole contribution is bigger than the continuum contribution. Releasing the above 
standard convergence criteria of OPE, we consider the ratio of perturbative contribution to the "total OPE 
contribution" (the sum of perturbative and other condensate contributions calculated) but not the ratio 
of perturbative contribution to each condensate contribution. Not bad, there are enumerably important 
condensate contributions (four-quark condensate and two-quark multiply mixed condensate) and other 
condensate contributions are much smaller than the perturbative contribution. Two important condensate 
contributions could cancel with each other to some extent, which brings that the ratio of perturbative 
contribution to the "total OPE contribution" is 71% at M 2 = 0.7 GeV 2 and increases with the M 2 . In 
this sense, the OPE converges when M 2 > 0.7 GeV 2 (note that the perturbative contribution of OPE 
series here is not always bigger than other terms in succession). Thus, the range of M 2 for K*K* is 
taken as M 2 = 0.7 ~ 1.3 GeV 2 for */io = 2.4 GeV. Similarly, the proper range of M 2 is obtained as 
0.7 - 1.4 GeV 2 for = 2.5 GeV, and the range of M 2 is 0.7 - 1.5 GeV 2 for = 2.6 GeV. In the 
chosen region, the corresponding Borel curve to determine the mass of K*K* is shown in the left panel 
of FIG. 2, and we extract the mass value 2.02 ± 0.11 GeV. In the end, we vary quark masses as well as 
condensates and arrive at 2.02 ± 0.11 ± 0.04 GeV for K*K* (the former error reflects the uncertainty due 
to the variation of so and M 2 , and the latter error is resulted from the variation of QCD parameters) or 
2.02 ± 0.15 GeV in a concise form. For KK* , we choose the minimum value of M 2 to be 0.7 GeV 2 in view 
of its OPE convergence. Furthermore, the ratio of pole contribution to continuum contribution from sum 
rule ([5]) for ^/sq = 2.2 GeV is approximate to 52% at M 2 = 1.3 GeV 2 . Thus, the maximum value of M 2 is 
taken as 1.3 GeV 2 for i/so = 2.2 GeV. With the similar analysis, the maximum M 2 is taken as 1.4 GeV 2 
for yfso — 2.3 GeV; for ^/lo = 2.4 GeV, the maximum M 2 is taken as 1.5 GeV 2 . The dependence on 
M 2 for the mass of KK* from sum rule ([6]) is shown in the right panel of FIG. 2. Finally, we arrive at 
1.85 ± 0.09 ± 0.05 GeV for KK* (the first error reflects the uncertainty due to the variation of sq and M 2 , 
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and the second error is resulted from the variation of QCD parameters) or 1.85 ± 0.14 GeV concisely. 

In Ref. it has been suggested that the two states Z~^ x and Z^ 2 in question should appear near their 
respective two meson thresholds, namely KK* and K*K*. We would make a comparison between the 
QCD Sum Rules' results here and the known thresholds. For the state, the result of QCD sum rule 
calculation here is approximately 320 ~ 600 MeV higher than the KK* threshold. For the Zf 2 state, our 
result is roughly 90 ~ 390 MeV higher than the K*K* threshold. 
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FIG. 1: In the left panel, the solid line shows the relative pole contribution (the pole contribution divided by 
the total, pole plus continuum contribution) and the dashed line shows the relative continuum contribution from 
sum rule (1101) for yfso = 2.4 GeV for K* K* . The OPE convergence is shown by comparing the perturbative, 
two-quark condensate, four-quark condensate, mixed condensate, two-quark multiply two-gluon condensate, two- 
quark multiply mixed condensate, six-quark condensate, mixed multiply mixed condensate, and four-quark multiply 
two-gluon condensate contributions from sum rule (|10[l for ^fs^ = 2.4 GeV for K* K* in the right panel. 
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FIG. 2: In the left panel, the dependence on M 2 for the mass of K* K* from sum rule (|11[) is shown. The 
continuum thresholds are taken as ^/so = 2.4 ~ 2.6 GeV. For y/so = 2.4 GeV, the range of M 2 is 0.7 ~ 1.3 GeV 2 ; 
for yfsE = 2.5 GeV, the range of M 2 is 0.7 ~ 1.4 GeV 2 ; for yfsE = 2.6 GeV, the range of M 2 is 0.7 ~ 1.5 GeV 2 . The 
dependence on M 2 for the mass of KK* from sum rule (|6]) is shown in the right panel. The continuum thresholds 
are taken as yfs^ = 2.2 ~ 2.5 GeV. For y/sE = 2.2 GeV, the range of M 2 is 0.7 ~ 1.3 GeV 2 ; for </sE = 2.3 GeV, 
the range of M 2 is 0.7 ~ 1.4 GeV 2 ; for ^/sE = 2.4 GeV, the range of M 2 is 0.7 ~ 1.5 GeV 2 . 



IV. SUMMARY AND OUTLOOK 



In theory, there could exist two charged strangeonium-like molecular states Z^ and Z^ 2 from an effective 
Lagrangian study. In this work, we have employed the QCD sum rule method to predict masses of Z^ and 
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Z^2, taking into account contributions of operators up to dimension ten in the OPE. Our final numerical 
results are 1.85 ± 0.14 GeV for Z+ (KK*) and 2.02 ± 0.15 GeV for Z+ (K*K*), which are both above 
their respective two meson thresholds. One can expect that these results could be helpful for investigating 
Z^ and Z^ 2 experimentally. We suggest to start the search for these two states in some decay process 
such as y(2175) — > 0(1O2O)7t + 7t~ in future experiments, especially Super-B, Belle II and BESIII. 
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